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ABSTRACT 

We have used archival FUSE data to complete a survey of interstellar HD in 
41 lines of sight with a wide range of extinctions. This follow up to an earlier 
survey was made to further assess the utility of HD as a cosmological probe; to 
analyze the HD formation process; and to see what trends with other interstellar 
properties were present in the data. We employed the curve-of-growth method, 
supported by line profile fitting, to derive accurate column densities of HD. We 
find that the N(HD)/2N(H 2 ) ratio is substantially lower than the atomic D/H 
ratio and conclude that the molecular ratio has no bearing on cosmology, because 
local processes are responsible for the formation of HD. Based on correlations with 
E(B-V), H 2 , CO, and iron depletion, we find that HD is formed in the densest 
portion of the clouds; the slope of the logN(HD)/log(H 2 ) correlation is greater 
than 1.0, caused by the destruction rate of HD declining more slowly than that 
of H 2 ; and, as a sidelight, that the depletions are density dependent. 
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1. Introduction 

The H 2 isotopologue HD was first detected in Copernicus spectra by Spitzer et al. 
(1973, 1974) and Morton (1975), and has subsequently been observed in many sightlines 
by FUSE. The HD lines arising from the lowest-lying rotational levels (J=0 and J=l) are 
far weaker than their counterparts for H 2 , but are in many cases still strong enough to be 
saturated, requiring a curve-of-growth analysis for column density deter minations. To date 



HD analyses from F USE spectra have been published only for a few stars (IFerlet et al. 2000 



Lacour et al. 20051 ) though many more detections reside in the FUSE archives. In addition, 
Lacour et al. have re-analyzed Copernicus spectra for several stars, resulting in a uniform 
survey of HD abundances in some 17 sightlines. The Copernicus stars generally have Ay 
values of 1.0 or less, while our FUSE targets cover a range of Ay from 0.1 to about 3 
magnitudes. 

The ratio of HD to H 2 — more specifically, N(HD)/2N(H 2 ) - in the Lacour et al. 
survey ranges from a few times 10~ 7 to several times 10 -6 . This is somewhat higher than 
the values found earlier from Copernicus data, but those values were based on two strong 
lines that were almost always saturated. 

We follow the convention of Lacour et al. (2005), and assert that, in a region where 
all hydrogen (including deuterium) is in molecular form, the atomic D/H ratio is the same 
as the N(HD)/2N(H 2 ). The reasoning is that the molecular fraction (N(HD)/N(H to t ai ) = 
N(HD)/[N(H I) + 2N(H 2 )]) reduces to N(HD)/2N(H 2 ) when H I can be neglected. For every 
HD molecule, there is one D atom, so the atomic ratio D/H is equal to N(HD)/2N(H 2 ). Note 
that this assumes that not only H I, but also atomic D, can be neglected. This simplifying 
assumption is justified in cloud cores, but it may break down in regions where hydrogen (and 
deuterium) is only partially in molecular form, which may include some of the less-reddened 
lines of sight in our survey. But even there, as shown in the next few paragraphs, we can 
not use HD as an indicator of the D/H ratio anyway. 

It would be useful if the N(HD)/2N(H 2 ) ratio were a reliable indicator of the atomic ratio 
D/H, because the latter is an important probe of the early expansion rate of the universe and 
therefore valuable in cosmology. An excellent summary of the observed atomic D/H ratio 
and its cosmological implications can be found in Linsky et al. (2006). But for three distinct 
reasons the comparison of HD to H 2 is difficult to apply to the cosmological problem. 

First, H 2 begins to be self-shielding much before HD, as a function of depth in a cloud. 
For example, H 2 becomes self- shielding at about N(H 2 ) = 10 19 cm~ 2 (which corresponds to 
N(HD) = 10 12-13 cm -2 ) and a similar column density would be required for HD - which is 
incredibly unlikely. For HD to become self-shielding, the column density of H 2 would have 
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to be of order 10 25 cm -2 ! Even infrared observations of the quadrupole transitions at 2.4 
jum would never reach this column density. The fact that HD is not protected from radiative 
dissociation acts to lower the N(HD)/2N(H 2 ) with respect to the atomic D/H ratio in the 
denser portions of the lines of sight. 

Second, if formation of HD occurs on grain surfaces (in parallel with H2 formation), the 
lower mobility of D atoms on grains as compared to H atoms, greatly reduces the formation 
rate of HD as compared to H 2 , acting to reduce N(HD)/2N(H 2 ) ratio relative to the D/H 
ratio. If grain formation of HD dominates, this would act to decrease the N(HD)/2N(H 2 ) 
ratio over the atomic ratio. For a more complex and complete analysis of HD formation on 
grain surfaces, see Cazaux et al. (2008). 

But third, HD has a gas-ph ase formation channel, through the " chemical fractionation" 



reaction H 2 + D + — > HD + H + fjWatson 1973] ; [Le Petit et al. 20021 ). w hich tends to enhance 



the N(HD)/2N(H 2 ) ratio and depends on the cosmic-ray ionization rate (IBlack fc Dalgarno 19731 ) 
The comparison of the N(HD)/2N(H 2 ) ratio to the atomic D/H ratio can tell us which pro- 
cesses are most important. 

In any case, the N(HD)/2N(H 2 ) ratio should not be expected to reflect the D/H ratio 
(Lacour et al. 2005). In general the N(HD)/2N(H 2 ) ratio is not useful for determining the 
cosmological D/H ratio, although it may be useful in constraining the rate of cosmic-ray 
ionization (which produces D + ). 

This paper, then, presents a new and more complete survey of HD than was avail- 
able previously, to test formation theories and possible correlations with other interstellar 
quantities. We have organized the paper into sections as follows: Observation and Data 
Reduction, Data Analysis, Determination of HD Column Densities, Correlations, Summary 
and Conclusions. 



2. Observation and Data Reduction 

The FUSE mission covered a wavelength range (905-1188 A) rich in ground state tran- 
sitions, making it one of the most useful tools to study the interstellar medium. The 41 
sightlines for this survey are archival FUSE spectra, chosen to span different cloud types and 
physical conditions. The sightline parameters used to select targets included color excess 
(0.11 < E(B - V) < 0.83), the indicator of grain size R„ (2.25 < R v < 4.76), molecular 
fraction (0.02 < f H2 < 0.76), and H 2 column density (19.07 < log N(H 2 ) < 21.11). Tabled 
lists all of the relevant physical properties for the sightlines included in this study. Another 
selection criterion was a high signal-to-noise ratio around the 1105.83 A HD line. This line 
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is 100 times weaker than the rest of the HD absorption features measured, all of which have 
similar oscillator strengths (HD f- values and wavelength data from Abgrall & Roueff 2006). 
In many cases this weak line's equivalent width is the only one falling on the linear portion 
of the curve of growth, thus providing unambiguous information about the column density 
as Lacour et al. found in their study. 



Because HD is known to be correlated with H 2 abundances (lLacour et al. 20051 ). our 
target list includes sightlines with H 2 column densities in the literature. H 2 column densities 
are taken from Shull et al. (in preparation, 2008) and Rachford et al. (2002; in preparation, 
2008). For this study the FUSE data (Table [2]) were pre-calibrated with version 3.1.4 or 
newer of the CALFUSE pipeline. For each observation we used a cross correlation analysis 
to align individual exposures with strong absorption features before co-adding them. 

In most cases we co-added the detector segments to increase the signal to noise ratio 
before measuring the HD equivalent widths. When combining the detector segments one 
significant systematic flaw was found in the data. The raw data from the LiF2A detector 
segment revealed the same systematic detector defect due to a Type I dead zone (as described 
in section 9.1.6 of The FUSE Instrument and Data Handbook) in all sightlines in the vicinity 
of the weak 1105.83 A HD line; thus we excluded that segment and only used LiFlB. 



3. Data Analysis 

To determine the column density of HD we measured equivalent widths for all possible 
HD lines (see Table H]) and performed a curve of growth (COG) analysis. All of the HD 
lines measured are transitions from the ground rotational state (J=0). While there are 
over 25 HD absorption lines in the FUSE range, only seven are isolated enough from other 
features to determine accurate equivalent widths. To measure the equivalent width we first 
defined a continuum on both sides of the feature and normalized it with low order (1 to 
3) Legendre polynomials. When there was no other absorption feature in the immediate 
vicinity we integrated over the HD line while summing the flux errors provided by FUSE in 
quadrature. If there was another absorption feature blended with the HD line, we fit the HD 
and interfering line with Gaussians (see next section for details). The errors in the width and 
depth of the Gaussian were propagated into the equivalent width error. The statistical error 
on the equivalent width was summed in quadrature with the continuum placement error. 
The continuum error was up to an order of magnitude smaller than the statistical error in 
all cases. 

Of the seven lines, six have similar oscillator strengths (see Tabled]); the seventh and 
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weakest line is located at 1105.86 A and is very important for constraining the COG fit. 
The weak HD line is difficult to measure in most cases and in other cases it is impossible 
to measure at all. Without a line to measure we used the signal-to-noise ratio to set a two 
sigma upper limit on the equivalent width using the equation: 

AAM 1 / 2 . . 

a = Js/W (1) 

where A A is the pixel scale, S/N is the signal-to- noise ratio, M is the width of the feature 
(here FUSE resolution elements (~9 pixels) are used since the width of the line will be 
smeared out by the resolution), and a gives the one sigma error of the equivalent width 
( 1 Jenkins et al. 19731 ). 



3.1. Modeling the C I* Line 

The first obstacle in measuring the weak HD line is a very crowded patch of continuum. 
In a study of the sightline to ( Oph, Morton (1978) found an unidentified absorption feature 
at 1105.92 A. In the current study we have been unable to confirm or deny the existence of 
a weak feature near this wavelength; thus, if such a feature does exist it could be a source of 
systematic error. Morton also rejected a feature at 1105.82 as being the weak line of HD due 
to its large equivalent width. The conclusion of Lacour et al. (2005) and the present study 
is that the feature at 1105.83 is indeed HD and its strength is consistent with all of the other 
measured HD features (see Figures [2J 0, [5] & M containing HD COG's). Shortward of the HD 
line, by 0.13 A, is an absorption feature from the first excited state of carbon as discussed 
by Lacour et al. In sightlines where the C I* and HD lines were blended we modeled the C 
I* and divided it from the spectrum before measuring the HD line. 

Removing this blended C I* line turned out to be more challenging than originally 
thought. The f- value of the 1105.73 C I* line, previously measured as 0.0113 by Morton 
(1978), did not fit the curve of growth set by the other C I* lines. This is probably due 
to uncertainty in the adopted f-value of this weak line. The previous value was based on a 
single measurement of the line in the spectrum of £ Oph. So to accurately model the C I* 
line we used existing archival data to make an empirical estimate of the f-value, based on 
more sightlines. 

To accomplish this task we found Copernicus, FUSE, and/or STIS data in sightlines 
where the 1105 A C I* line was well defined, not heavily blended with other lines in the 
region (HD 12323, o Per, HD 207538, HD 210839) and measured the equivalent widths of 
all well defined C I* lines, excluding the 1105 A line. For each sightline we then fit these 
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equivalent widths to a single component curve of growth varying both 6-value and column 
density (see Figured] for an example). Once a best fit curve of growth was found, the f- value 
of the 1105 A line was calculated using its measured equivalent width. Uncertainty in the 
f- value was assigned based on errors in the measured equivalent width of the 1105 A line. 
Calculated f- values and uncertainties for each sightline can be found in Table [61 The results 
from all sightlines were combined using a weighted average and the final f-value was found 
to be 0.00621° °°^. 

The C I* lines used in modeling are in Table [51 We used STIS data when available, 
along with FUSE data to measure the C I* absorption features. Using the best fit column 
density and 6-value of C I*, we generated a Voigt profile of the 1105.73 A C I* line. This 
profile was then convolved with the Gaussian instrumental profile (assuming a resolution of 
20 km s -1 ) and divided out of the spectrum, allowing us to measure the HD 1105.83 A line. 

We began the current study using C I* line f- values from Morton's compilation (1991). 
However, in light of newer theoretical f-values from Zatsarinny & Froese Fischer (2002) and 
Froese Fischer (2006), it was important to determine whether differences due to revised f- 
values would substantially affect our HD column densities. Repeating the calculation of the 
C I* 1105.73 A line f-value yielded only a very small change (about three percent) from 
the previously calculated f-value. All C I* column densities were also recalculated using the 
newer f-values. Total column densities did not vary greatly from those using the Morton 
(1991) f-values. In all but three cases differences in C I* column density using the two sets 
of f-values differed by 0.1 dex or less. The three most discrepant cases were HD 73882, HD 
101436 and HD 149404. For these three cases the equivalent width of the weak HD line 
was remeasured and found to be well within the la errors of the previous measurements. 
Similarly, in all three cases rerunning the COG fit yielded only small changes in the HD 
column density that were within the la error bars of the previous measurements (in all cases 
the difference was less than 0.1 dex). Thus, we did not think it necessary to repeat all of the C 
I* modeling, and remeasure the weak HD line, as other errors far exceed any error introduced 
by minor changes in some of the C I* f-values, and changes in C I* column densities were 
not substantial in the vast majority of cases. All reported C I* column densities are those 
using the newer f-values from Zatsarinny & Froese Fischer (2002) and Froese Fischer (2006). 



4. Determination of HD Column Densities 

Once the influence of C I* had been removed, we were able to do a curve-of-growth 
analysis to determine the column density, but we used profile fitting in a few Cases as 8b 
check. We did two different COG analyses; multiple component COG's for a sub-set of 13, 
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and single component COG's for all 41 sightlines. Both are described below in subsection 
4.1 'Curve of Growth Analysis'. The profile fitting is described in subsection 4.2. 



4.1. Curve of Growth Analysis 

In many cases multiple cloud components exist along a single line of sight increasing the 
difficulty in accurately measuring the column density. Having knowledge of such structure 
is most important when the equivalent widths fall on the flat portion of the curve. We can 
create a multiple component COG, but first we need to find a suitable tracer of HD. 



HD has been found to be correlated to H2 ( lLacour et al. 20051 ). Similarly CH is traced 
by H2 because of the chemistry: the creation of CH in the diffuse ISM is directed through 
a series of gas phase reactions with H 2 (Danks, Federman & Lambert 1984; van Dishoeck & 
Black 1989; Magnani et al. 1998). Since CH and HD are both correlated to H 2 abundances, 
there is an inferred correlation between CH and HD, which has been verified by Lacour et 
al. 

For thirteen sightlines we had high resolution CH data (Welty, private communication) 
which we used to define a velocity structure for HD. From the measured CH line-of-sight 
velocity structure, we generated multiple component COG's. The HD equivalent widths were 
fit to the curve varying only the total column density. For comparison the data were also 
fit to a single component COG that found a best fit of both 6-value and column density by 
minimizing the x 2 ■ Figures [2] and [3] show the side by side comparison of the COG analyses, 
and Figure H] compares the column densities and errors from the two methods. 

The COG's from the two methods fit the data in all but three cases. The side by side 
comparison show a clear mis-fit of the multiple component COG for HD 149404, HD 185418, 
HD 192639, even though the column densities are within one-sigma errors for the first two 
targets. The discrepancy arose because a multiple component COG assumes the best fit 
velocity structure of CH is the actual structure of the HD and does not allow for errors. 
Varying the 6-values on the multiple component COG within their one sigma errors gave 
solutions that were consistent with the data and the single component COG. 

With the reasonable assurance that our single component COG model is good we ac- 
counted for any other systematic errors that we could not quantify by setting the reduced 
X 2 of the best-fit curve equal to one. In 22 sightlines the reduced x 2 was already less than 
one, thus the errors were unchanged; for the remaining 19 cases the assumed good fit scaled 
up the errors in equivalent width and column density. 
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The consistency between the results of the single and multiple component COG demon- 
strates that a single component model is a good approximation when we do not have compo- 
nent information. Thus, we are justified in combining the HD results for the 13 targets with 
component information and the 28 targets without; see Figures and O for single component 
COG's and Table [7] for column densities. 

Solving for column density in the COG analysis was more complicated when we only 
had an upper limit on the 1105 A line. We used the limit on the weak line to constrain the 
column density by throwing out solutions beyond the upper limit. The limit was not used 
in calculating the \ 2 °f the fit. 

Lower limits on the b- value were determined when all the points fell on the linear portion 
of the COG. The 6-value controls where the curve will transition from the linear to the 
saturated portion. Since the transition cannot occur lower than the data points a two sigma 
lower limit is set (HD 91824, HD 93204, HD 93205, HD 94493, HD 161807, and HD 201345.) 

4.2. Profile Fitting 

Profile fitting of the HD lines was done for five lines of sight to verify the COG column 
density measurements. All the available HD lines were fit simultaneously so that a single best 
fit column density and b- value would be measured. First the continuum surrounding each HD 
feature was normalized by fitting a low order Legendre polynomial to the continuum on both 
sides. Shifts in the central velocity of the HD features were noticed from one line to another, 
(with each detector segment showing different velocity shifts, usually on the order of ~10 km 
s _1 ) so the lines were each fit with a single Gaussian and the central wavelength corrected so 
that each feature would be at its rest wavelength. Profiles were modeled assuming a single 
Voigt profile with column density N, line width b and central velocity v. These profiles were 
then convolved with a Gaussian instrumental profile. We used the non-linear least-squares 
curve fitting algorithm MPFIT by C. Markwardt0 to find the best fit parameter values. 
MPFIT is a set of routines that uses the Levenberg-Marquardt technique to minimize the 
square of deviations between data and a user-defined model. These routines are based upon 
the MINPACK-1 Fortran package by More' and collaborator^. 

Attempts were also made to fit multiple components when previous measurements of CH 
existed. However, even when the velocity structure was fixed by the CH structure, individual 



^ttp : / / cow . physics . wise . edu/~craigm/idl/idl . html 
^http : //www . netlib . org/minpack 
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component column densities did not converge to reasonable answers. Usually one component 
of the fit dominated while the other component (s) were an order of magnitude or two lower 
in column density (differing greatly from the relative strengths of the CH components). 
Because of this, we report only the single component fits in Table El 

For all five sightlines we find that the column densities and 6-values for both the curve 
of growth and profile-fitting method are within standard errors of one another (including 
HD 192639, our previously discrepant case). We also find that the best fit profiles match the 
absorption profiles in the data well (see Figure [7j). Thus, we are confident that our column 
densities measured through curve of growth fitting are reliable. 

4.3. Comparisons with previous results 

We are able to compare our column densities with the previous HD survey of Lacour et 
al. (2005). This survey had 7 stars observed with FUSE and 10 more lines of sight observed 
years before, with Copernicus. The comparison of our results are in Table [31 

We re-analyzed the 7 FUSE targets that Lacour measured to check whether our method 
was consistent with their results. Our method differs in a few small ways. We re-derived the 
C I* f- value (as described above), so our de-blending of the weak HD line is different. Also 
we weighted each line equally in the COG fit, while the Lacour group doubly-weighted the 
optically thin HD line. Even with these difference our column densities are consistent within 
the uncertainties. 

We also did single-component curves of growth for every sightline, which the Lacour 
group did not, and we found consistent column densities with the multicomponent results. 
This allowed us to expand our survey with confidence, in the end obtaining accurate column 
densities for 41 stars. 

Previous HD measurements made by Spitzer et al. (1974) are included in Figure [S] and 
the correlation analysis. The plot includes the upper limits for 5 Per and A Ori and the 
measurements for n Sco, 5 Sco, o Sco, and 7 Ara. 

5. Correlations 

One way to derive information about the formation mechanism of HD, and to test 
theoretical chemical models, is to see how various interstellar quantities correlate with HD 
column densities. Before we do that, a few words of caution are in order. Of necessity, we 
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can measure only integrated line-of-sight quantities, not localized to the same portion of the 
lines of sight. For example, any correlations of molecular quantities versus dust indicators 
or H I are probably not informative, because both dust and H I are everywhere, inhabiting 
every region from diffuse to dark clouds. We assume, with good reason, that the observable 
HD and H2 (along with CH) are confined to diffuse molecular and translucent clouds (as 
defined by Snow & McCall 2006), so we should not expect to find good correlations with 
either dust or H I. If we do find good or decent correlations, that would indicate that diffuse 
molecular and translucent clouds dominated the lines of sight in our survey. On the other 
hand, those quantities that do correlate well with each other would show that they arise in 
the same portion of the line of sight. All correlations can be found in Table [9l including 
some not illustrated. 

The first correlation of interest is the one between HD and H 2 , which is shown in Figure 
[HJ There we see a very good correlation, suggesting that the formation of HD depends on the 
local abundance of H2, or that both are formed by the same or similar processes. Because H2 
shoots up in column density due to the self-shielding about when N(H2) reaches 10 19 cm" 2 , 
we checked for a change in slope around that point. We see more scatter below log N(H2) 
= 19.5 and no significant change in slope (below 19.5, the slope is 1.55, based on a fairly 
narrow range; above, it is 1.33). Considering all of our data points together gives an overall 
slope of 1.25 and a correlation coefficient is 0.94. 

The significant departure from a slope of 1.00 in this correlation can be interpreted in 
two different ways: either HD is being formed at a faster rate than the formation rate of H 2 in 
the portion of the clouds where both are forming; or HD is destroyed at a lesser rate than the 
destruction rate of H2 in those regions. This latter explanation may be preferred, because 
of optical depth effects and how they affect the photodestruction rates. For most of the 
observed range in Ay, the H 2 J = and 1 lines are damped, while the corresponding lines of 
HD lines are becoming saturated - but not damped. In this range in Ay (between Ay about 
0.1 and 0.8), the rate of HD des truction is rapidly de clining, while the photodestruction rate 



of H2 is only slightly declining (ILe Petit et al. 20021 ). This situation mimics the appearance 



of relative rate of formation favoring HD, and that explains the slope of the N(HD)/N(H 2 ). 

The alternative explanation of the slope in Figure [HI that HD is being formed faster 
than H2, could only be viable if some additional formation process were at work, in addition 
to the standard gas-phase reaction network. And that would imply HD formation on grain 
surfaces competes with the gas-phase mechanism. But the detailed analysis of HD formation 
on grains by Cazaux et al. (2008) explores conditions in which grain formation of HD could 
be important, and our observed lines of sight do not meet the criteria outlined by Cazaux et 
al. So our explanation, that HD has started to be destroyed at a lesser rate, is more likely. 



- 11 - 



This finding of a greater slope than 1.00 seems to be inconsistent with the gas-phase 
models of Le Petit et al. (2002), who predict that the ratio of HD over H 2 should decline over 
our entire range in molecular fraction. Their Figure 4 shows the ratio N(HD)/2N(H 2 ) as a 
function of the molecular fraction of hydrogen, and this figure shows a steady decline with 
increasing molecular fraction. Only above a molecular fraction of 0.9 does the ratio start to 
increase. The study by Le Petit et al. only considers a single cloud, whereas real sightlines 
cover many clouds in most cases. But it is difficult to see how this helps reconcile the model 
with the observations, because the situation only gets worse when multiple clouds, which 
collectively form the total column density, are considered. The model would be correct only 
if all the molecular regions along the lines of sight, each having a molecular fraction > 0.9, 
contained all the HD and H2, with no spillover to less dense regions. 

Figure shows the correlation between the column densities of HD and CO (the CO 
data were taken from Burgh et al. 2007). Here we see a very good correlation (coefficient 
0.93), which we take as an indication that related processes account for both HD and CO. 
Since CO is formed by gas-phase reactions (e.g Kaczmarczyk 2000a; 2000b; and references 
cited therein), we again conclude that HD is also likely formed by gas-phase reactions. The 



ion-neutral reaction H 2 + D + — > HD + H + (jLe Petit et al. 20021 ). with the ion D + being 



created by cosmic-ray impacts, is commonly thought to be the primary gas-phase formation 
channel for HD. 

To test our assertion that HD and H 2 do not inhabit the whole line of sight in most 
cases, we show the correlations of both with quantities that are expected to arise over the 
entire sightline. In Figure [TUl we show the correlations of HD with the dust parameters total 
H (i.e., N(H I) + 2N(H 2 )); Figure [H] shows the correlations of H 2 and HD with E(B-V); and 
Figure [T2l shows the correlations of H 2 and HD with Ay. These correlations do not differ 
significantly from each other, except there is a hint that both HD and H 2 go with E(B-V) 
slightly better than they correlate with either total H, or Ay. Neither HD or H 2 correlate 
significantly with Ry. 

The correlations of HD and H 2 with E(B-V) are better than those with H I, telling 
us something surprising and perhaps useful: that the grains responsible for differential ex- 
tinction tend to be in the densest parts of the observed lines of sight, where HD and H 2 
are present. The grains responsible for E(B-V) are thought to be comparable in size to the 
wavelengths being scattered, so we conclude that other grains in the lines of sight may have 
a different size distribution. 

We tested correlations of both H 2 and HD against the molecular fraction, given by f 
= N(H 2 )/N(Hr ota ;), and found a fairly good correlation (see Table [9]). Neither should be a 
surprise, because they just show that H 2 and HD exist primarily in the densest portions of 
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the lines of sight, as expected. 

Finally, we have looked at correlations between iron depletion, as shown by the ratio of 
iron to total hydrogen, and we see several interesting things. These correlations are shown 
in Figure [13j We could have considered correlations with depletions of other elements, but 
we have more complete data on iron for our lines of sight, and we expect other refractory 
elements to follow more or less the same trends as iron. The iron column densities come 
from Jensen & Snow (2007), and are listed in Tabled) 

The three correlations we considered are consistent with the idea that the depletions 
are density dependent. Remembering that an increasingly negative value of the Fe/H ratio 
means more depletion of iron, we see that the depletion is growing with indicators of density 
such as the HD column density, the fractional ratio of HD (i.e., N(HD)/N(H<r ota /); and the 
ratio of HD to H 2 (see Figure [TBI) . The indicators we see here, showing that depletion 
increases with density, is consistent with the finding of Burgh et al. (2007), who found that 
the depletions correlate with CO, another density indicator. 



6. Summary and Conclusions 

We have derived column densities of interstellar HD for 41 lines of sight, covering a wide 
range of total extinctions, total hydrogen column density, and differential extinction E(B-V). 
We went to considerable effort to reduce error in the derived column densities, which were 
ultimately found using the curve-of-growth method, after measuring the equivalent widths. 
One challenge we faced, and successfully resolved, was that the weakest line of HD, detectable 
only toward highly-reddened stars, was contaminated by a weak C I* line. We had to model 
and remove the C I* line in those cases. Where possible, we checked our results against 
other published column densities, with generally good agreement. 

We found values of the logarithmic ratio of N(HD)/2N(H 2 ) ranging from -6.18 to -5.13 
- generally much lower than the galactic atomic D/H ratio and even lower with respect to 



the extragalactic (or the primeval) ratios (ILinsky et al. 20061 ) . This says two things: (1) the 



enhanced photodissociation of HD over H 2 dominates and lowers the N(HD)/2N(H 2 ) ratio; 
and (2) our initial assumption that we can not use N(HD)/2N(H 2 ) ratio as a cosmological 
test is confirmed, and only the atomic ratio gives a true value of the primordial production 
of deuterium (if that; see Linsky et al. 2006). 

In order to see whether our results are consistent with extant chemical models, we 
examined correlations of HD with other interstellar parameters, finding the following: (1) 
HD is declining more slowly than the destruction rate of H 2 ; (2) N(HD) correlates well with 
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N(CO); and (3) the depletion of iron (and probably other refractory elements) is enhanced 
in the densest portions of the lines of sight. The model of Le Petit et al. (2005) of the 
formation of HD is inconsisten with the observations, unless individual dense clouds along 
the lines of sight contain all the HD and H 2 . 

We are very grateful to the anonymous referee, who stimulated us to revise this paper 
in many useful ways. We are grateful to several colleagues (including especially Dan Welty 
and Mike Shull) who gave us their advice, not to mention their data; and to NASA, which 
funded this research through grants NNG05GA85G and NAS5-98043. Also thanks to Ms. 
Lynsi Aldridge who helped with the early analysis of the data. 
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Table 1. Sightline Properties 



Star 


E(B-V) 


Rv 


Ref. a 


Molec. 


log N(H 2 ) 


log N(CH) 


log N(H Toi ) 


Rcf. a 


log N(Fe II) b 


log N(C I*) c 


(HD) 








Fract. 


(cm" 2 ) 


(cm" 2 ) 


(cm" 2 ) 




(cm" 2 ) 


(cm" 2 ) 


12323 


0.26 


2.75 


1 


0.21 
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A 
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14. yb ± U.UZ 


Z777o 
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3 
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45314 


0.46 


4.42 


1 


0.42 


20.60 




21.28 


4 






53367 


0.74 


2.38 


3 


0.52 


21.06 




21.65 


3 


15.08 


15.1o ± l.bo 


(oool 


0.72 


3.37 


3 


0.67 


21.11 


13.57 


21.58 


5 


15.15 


lb.o5 ± U.b4 


( 4yzu 


U.oo 






n 01 
U.Z1 


on ofi 
ZU.Zb 
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Z1.Z5 


A 

4 






nnn07 

yuuo r 


0.30 






0.08 


19.77 




21.18 


4 


15.15 




yiooi 


n qo 
U.oU 


0.05 


i 

i 


n no 
U.UZ 


1 n H7 

iy.u r 




01 1 P. 


A 

L 


1 K OQ 

15. Zo 




91824 


0.26 


3.35 


1 


0.09 


19.84 




21.19 


4 






noon a 
yozU4 
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0.04 
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A 
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O.Z5 


i 
i 


n ne; 
U.U5 
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iy. / o 
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A 

L 


15. oo 
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14. lo ± U.U5 


yozUb 


0.39 






0.03 


19.52 




21.35 


A 

4 




1 A 1 Q _l_ n 1 *7 

14. lo ± U.17 




0.36 


4.7b 


1 


0.03 


1 n 11 

iy. 7 ( 




21.55 


A 

4 




i q l 7 /r 7 _L n n /i 
lo. / ( ± U.U4 


94493 


0.20 






0.18 


20. 12 




21.17 


A 

4 


15.38 




101131 


0.34 








20.27 






A 

4 






1 ni 1 on 


U.oo 


O A Q 

Z.4o 


i 


n o7 
U.Z ( 


on a o 
ZU.4Z 




01 on 
zi.zy 


A 

L 




i /i c/i i n 1 1 
14.54 ± U.ll 


101413 


0.36 






0.22 


20.38 




21.34 


4 






i ni a q« 
1U1400 


U.oo 






n oo 
u.zz 


on qo 
ZU.oo 




01 Q A 


4 




1 A 4f\ 1 n oo 

14. 4y ± U.ZZ 


1U4 1 U5 


n ofi 

U.ZO 


O Q1 

Z.bl 


i 


n i o 
U. i-Z 


1 n tiQ 

iy.yo 




oi i a 
Zl.lo 


A 

L 


1 r oo 
15. ZZ 
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n a n 
U.4U 






U.OO 


on fi/i 
ZU. 04 


iq in 

lo. iy 


01 on 
Zl.ZU 


5 




i/i t; c _l n ofi 
14.05 ± U.ZO 


lluooZ 


u.zz 


O /I o 
Z.4Z 


i 


n i o 
U. i-Z 


1 n 7Q 
iy. lo 




01 m 
Z1.U1 


/J 
4 






1 /I 7Q Q Q 

14 (boo 


0.52 


4.06 


3 


0.10 


20.46 




21.71 


3 


14.89 


i c no _L n i n 
15. Uz ± U.1U 


1 A Q A OO 

1484ZZ 


n oq 
U.Zb 


q no 
o.UZ 


i 

i 


n i a 
U. ID 


on 1 q 
ZU. lo 




01 OQ 

Zl.Zo 








1 /i O/i ri/i 
14y4U4 


n fiQ 

U.OO 


Q OQ 

o.Zo 


Q 
O 


n qq 
U.oo 


on 7n 

zu. /y 




01 c;7 
Z1.5 ( 


Q 

o 


1 OQ 

lo. Zo 


i a ok _i_ n no 

14. ob ± u.uy 


152233 


0.45 


2.95 


1 


0.17 


20.29 




21.37 


4 






152248 


0.45 


3.68 


1 




20.29 






4 






152723 


0.42 


3.36 


1 


0.13 


20.29 




21.49 


4 






161807 


0.11 








19.86 






4 






177989 


0.25 


2.83 


1 


0.23 


20.12 




21.06 


4 


14.81 


14.52 ± 0.03 


185418 


0.51 


2.32 


3 


0.47 


20.76 


13.12 


21.56 


5 




14.35 ± 0.08 


192639 


0.66 


2.84 


3 


0.32 


20.69 


13.45 


21.47 


5 




14.54 ± 0.08 


199579 


0.36 


2.95 


3 


0.38 


20.53 


13.36 


21.25 


5 




14.78 ± 0.12 


201345 


0.17 






0.03 


19.43 




20.91 


6 






206267 


0.53 


2.67 


3 


0.42 


20.86 


13.41 


21.54 


5 




14.91 ± 2.09 


207198 


0.62 


2.42 


3 


0.38 


20.83 


13.56 


21.55 


5 




15.18 ± 0.16 


207538 


0.64 


2.25 


3 


0.43 


20.91 


13.63 


21.58 


5 




15.23 ± 0.15 


224151 


0.42 


2.64 


2 


0.26 


20.57 




21.45 


4 


15.38 


14.34 ± 0.06 


303308 


0.43 


3.02 


1 


0.11 


20.24 




21.50 


4 


15.46 


14.38 ± 2.61 


308813 


0.28 






0.22 


20.29 




21.26 


4 







a References: (1) Valencic et al. 2004; (2) Rachford, private communication; (3) Rachford in preparation, 2008; 
(4) Shull in preparation, 2008; (5) Rachford ct al. 2002; (6) Burgh et al. 2007 

b Jensen & Snow 2007 
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Table 2. Observation data 



Star 


Stellar 


V 


Dist a 


Rcf. 


FUSE Data ID 


STIS Data ID 


(HD) 


Type 


(mag) 


(PC) 








12323 


09V 


8.92 


3900 


1 


P1020202000 




15558 


05IIIf 


7.95 


2187 


2 


P1170101000 




24534 


09.5pe 


6.10 


398 


2 


P1930201000 


o66p01020 


27778 


B3V 


6.33 


220 


3 


P1160301000 


o59s01010 


45314 


09pe 


6.60 


799 


2 


P1021301000 




53367 


BoIV:e 


6.94 


247 


3 


P1161101000 




73882 


09III 


7.24 


925 


1 


P1161302000 




74920 


08 


7.54 


1497 


2 


P1022601000 




90087 


09.5111 


7.80 


2716 


2 


P1022901000 




91651 


09Vn 


8.87 


3500 


1 


P1023102000 




91824 


07V 


8.15 


4000 


3 


A1180802000 




93204 


05Vf 


8.48 


2630 


2 


P1023501000 




93205 


03V 


7.76 


2600 


1 


P1023601000 


o4qx01030 


93206 


09.7Ibn 


6.24 


2512 


2 


P1023401000 




93222 


07IIIf 


8.11 


2900 


1 


P1023701000 


o4qx02030 


94493 


B0.5Iab 


7.27 


3327 


2 


P1024101000 




101131 


06Vf 


7.16 


709 


3 


P1024901000 




101190 


06Vf 


7.27 


2399 


2 


P1025001000 




101413 


08V 


8.33 


2399 


2 


P1025301000 




101436 


06.5V 


7.58 


2399 


2 


P1025401000 


o61z51010 


104705 


BOIb 


7.83 


3898 


2 


P1025701000 




110432 


BHIIe 


5.24 


301 


3 


P1161401000 




116852 


09III 


8.49 


4760 


2 


P1013801000 




147888 


B3 


6.78 


136 


3 


P1161501000 


o59s05010 


148422 


Blla 


8.64 


8836 


2 


P1015001001 
P1015002001 
P1015003001 




149404 


09Iae 


5.47 


1380 


2 


P1161702000 




152233 


06III 


6.59 


1905 


2 


P1026702000 




152248 


07Ibnfp 


6.11 


1758 


4 


P1026801000 




152723 


06.5IIIf 


7.31 


1905 


2 


P1027102000 




161807 


BOIIIn 


7.01 


383 


3 


P1222302000 




177989 


B0III 


9.34 


4909 


2 


P1017101000 




185418 


B0.5V 


7.45 


950 


5 


P1162301000 


o5c01q010 


192639 


07Ib(f) 


7.11 


1100 


1 


P1162401000 


o5c08t010 


199579 


06Ve 


5.96 


794 


2 


P1162501000 




201345 


09V 


7.75 


1907 


2 


P1223001000 




206267 


06.5V 


5.62 


850 


1 


P1162701000 




207198 


09.5Ib-II 


5.96 


832 


2 


P1162801000 


o59s06010 


207538 


09V 


7.31 


832 


2 


P1162902000 


o63y01010 


224151 


B0.5II 


6.05 


1355 


2 


P1224103000 


O54308010 


303308 


03Vf 


8.21 


2630 


2 


P1221602000 


o4qx04010 


308813 


09.5V 


9.32 


2398 


2 


P1221903000 





a Reference: (1) Savage et al. 1985; (2) Diplas & Savage 1994; (3) The Hipparcos 
and Tycho Catalogues; (4) Raboud ct al. 1997; (5) Sonncntrucker ct al. 2003 
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Table 3. Previously Reported HD Column Densities Comparison 



Previous Values 3, Our Values 
Star log N(HD) log N(HD) 

(HD) (cm- 2 ) (cm- 2 ) 



27778 15.51 to.33 15.84 ± 0.12 

73882 15.76 to.H 16 - 03 ± °- 34 

110432 15.28 to.u 15 - 30 ± °- 08 

185418 15.63 +° \l 15.51 ± 0.32 

192639 15.18 t°' 2 ° 15.57 ± 0.22 

206267 15.32 to ll 15.54 ± 0.10 

207538 15.70 t„ 21 15.82 ± 0.09 



a Valucs from Lacour et al. 2005 
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Table 4. Equivalent Widths of HD lines 



Star 


W959.82 


W975.58 


Wioil.46 


W1021.46 


W1054.29 


W1066.27 


W1105.86 


(HD) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 



1 9999 


43.9 ± 3.0 


53.0 ± 3.8 


60.4 ± 3.3 


55.1 ± 2.5 


^ n -i- 9 9 
00. u nz z.o 


46.4 ± 3.0 


7 J. 9 c 
0.1 1 z.o 


10000 




















86.0 ± 19.0 


on 1-1-19^ 
yu.i ztz iz.o 


70.3 ± 10.2 


9ft _i_ 7 « 
zo.y zn t.o 


9/1 ^94 
Z^tOO^i 


25.8 


± 5.2 


26.5 ± 10.7 


26.2 ± 4.2 


28.4 ± 2.8 


97 9 -l- i S 
z / .0 ztz 1.0 


32.3 ± 2.3 


10 j. 1 4 
lo.o in ±.ft 


9777S 
z it 1 














34.7 ± 5.6 








ZtZ Z. 1 


32.9 ± 4.4 


99 n J. 9 9 
zz.u ztz z.z 


/I KOI A 














57.0 ± 1.7 


59.9 ± 1.4 


en 9 J. 1 9 

oy.z zn. ±.z 


53.4 ± 1.1 


/ or; 7 
<_ ZO. t 


^99fi7 
OOOO f 


53.8 


± 10.5 














59.8 


± 15.8 


oi.y ztz 0. 1 


61.2 ± 8.6 


9n 4- 1 Q 
zo.o ztz ^.y 


79S&9 




















43.0 ± 10.1 


A A 9 J- fi Q 
<4:4:.Z ZtZ 0.0 


45.6 ± 16.0 


oa n -i- s n 

Z4:.U ztz o.u 


7/1090 

/ ^tyzu 


37.0 ± 3.3 


38.9 


± 


3.8 


43.7 ± 1.9 


42.2 ± 1.5 


49 4-1-17 
'-tz.^t zn 1.1 


36.0 ± 1.6 


11 n 4 n 8 
10. u ztz U.O 


00087 


20.6 ± 2.0 


26.1 


± 


3.6 


30.3 ± 1.4 


30.7 ± 1.5 


99 7-1-19 
ZO. 1 I LZ 


17.5 ± 1.1 


s~ A Q 

^.y 


Q1 fiM 

y 1001 


3.2 ± 


1.4 








5.3 ± 0.9 


6.2 ± 1.3 


9 q -I- n s 
o.y ztz u.o 


2.9 ± 1.0 




y loz^t 


15.4 


± 


2.3 


19.1 


± 


2.6 


34.1 


± 


6.9 


33.7 


± 


2.3 


91 7 _L 1 7 
Zl. I ztz 1.1 








<- O.O 


09904 


19.8 


± 


5.0 


23.9 


± 


4.1 


41.5 


± 


4.4 


32.4 


± 


6.8 


OQ 1 -LOO 
ZO. 1 Z\Z z .z 










0990^ 

yozuo 


19.8 


± 


1.3 


26.6 


± 


1.5 


37.6 


± 


1.4 








97 A -4- 1 n 
z ( .t± zn 1 .u 








1 ft _|_ fi 

1.0 ztz u.o 


Q990fi 


14.8 


± 


1.5 


18.3 


± 


2.2 


23.2 


± 


1.5 








17 7 4- fl 

1 l . 1 Z\Z U.O 


14.1 


± 


0.8 


/ Q 1 
<^ 0. 1 


09999 


14.8 


± 


2.1 


17.4 


± 


2.4 


25.4 


± 


1.9 


24.4 


± 


1.6 


90 1 _|_ 1 Q 
ZU.l it l.O 


13.3 


± 


1.8 


ft ft 
<. O.O 


Q/1/IQ9 


13.3 


± 


2.0 


16.2 


± 


3.6 


20.6 


± 


1.6 








17 9-1- 1 1 
1 I .0 ZIZ 1.1 








<- O.Z 


101191 
lUllol 


19.7 


± 


2.7 


24.5 


± 


2.9 


25.7 


± 


2.2 


26.6 


± 


1.7 


9R 9-1-19 

zo.o m 1.0 


24.7 


± 


1.5 


A 9 -4- 1 A 

*±.z zn. ±.*± 


1011 Qfl 

iui iyu 


22.1 


± 


2.6 


23.3 


± 


2.4 


27.1 


± 


2.0 


28.0 


± 


2.1 


9R 1-1-11 

zo. 1 zn 1.1 


26.6 


± 


1.4 


-I- 1 n 
o.u ztz l.U 


101413 


20.7 


± 


4.1 








25.4 


± 


3.7 


25.2 


± 


2.4 


9C 7 _|_ 9 K 

zo. / ztz z.o 


21.1 


± 


3.6 


ft _|_ 9 9 
0.0 ztz z.z 


i 01 /L9fi 
lUl^toO 


18.9 


± 


3.1 


20.2 


± 


2.9 


20.4 


± 


2.6 


21.0 


± 


2.8 


90 Q -4- 1 A 

zu.y in l.*± 


18.0 


± 


1.6 


m 9 4 1 n 

1U.Z ztz l.U 


104705 


9.4 ± 1.5 


12.6 


± 


2.9 


13.8 


± 


1.2 


14.0 


± 


1.6 


13.0 ± 0.7 








< 3.9 


110432 


25.9 ± 4.0 


27.1 


± 


6.2 


30.4 


± 


2.8 


30.2 


± 


2.1 


34.8 ± 2.1 


32.1 


± 


1.7 


13.2 ± 1.5 


116852 


8.9 ± 1.2 


13.0 


± 


1.4 


15.3 


± 


0.9 








13.7 ± 0.7 


10.6 


± 


0.8 


< 6.1 


147888 














15.6 


± 


4.8 


15.1 


± 


6.0 


15.6 ± 2.4 


14.2 


± 


3.0 


8.8 ± 1.6 


148422 


23.5 


± 


5.0 


21.9 


± 


8.5 


22.9 


± 


3.3 


26.2 


± 


4.6 


26.0 ± 2.6 








7.3 ± 2.1 


149404 


47.3 


± 


7.0 


56.7 ± 11.1 


61.1 


± 


7.4 


53.7 


± 


3.5 


63.7 ± 3.9 


55.5 


± 


4.6 


26.4 ± 2.5 


152233 


28.5 


± 


1.9 








31.1 


± 


1.8 


35.4 


± 


1.7 


32.0 ± 1.5 


29.1 


± 


2.1 


6.5 ± 1.6 


152248 


27.8 


± 


3.9 








33.3 


± 


2.2 


35.8 


± 


1.8 


30.6 ± 1.8 


28.2 


± 


2.2 


9.3 ± 2.1 


152723 








41.2 ± 3.5 


46.1 


± 


6.1 


46.3 


± 


3.1 


49.5 ± 2.4 


43.0 


± 


1.9 


6.6 ± 1.4 


161807 


13.5 


± 


3.0 


18.3 ± 4.6 


32.1 


± 


2.7 


27.7 


± 


2.7 


17.4 ± 1.8 


15.3 


± 


1.8 


< 6.9 


177989 


25.5 


± 


3.6 


33.3 ± 


8.8 


33.6 


± 


2.6 








26.5 ± 1.9 


26.6 


± 


2.7 


6.3 ± 1.6 


185418 


32.0 


± 


6.8 








46.3 


± 


6.4 








50.8 ±3.9 


49.7 


± 


4.2 


17.4 ± 4.9 


192639 


48.2 


± 


7.8 








51.9 


± 


6.1 


54.2 


± 


4.0 


54.7 ± 3.7 


53.4 


± 


4.1 


20.1 ± 4.0 


199579 


52.6 


± 


4.3 


52.7 ± 5.5 


71.6 


± 


4.5 


64.6 


± 


5.3 


61.9 ± 2.6 


60.3 


± 


2.7 


8.4 ± 1.2 


201345 


11.5 


± 


0.8 


15.9 ± 1.4 


22.2 


± 


1.0 


22.3 


± 


1.0 


14.8 ± 0.8 


10.6 


± 


0.8 


< 2.4 


206267 


42.5 


± 


2.5 


48.5 ± 9.1 


46.1 


± 


5.1 


52.0 


± 


3.7 


49.5 ± 2.4 


45.6 


± 


3.2 


24.1 ± 2.1 


207198 


47.1 


± 


7.8 


53.8 


± 12.3 


57.8 


± 


7.1 


47.3 


± 


4.6 


56.9 ± 2.8 


51.1 


± 


4.3 


34.7 ± 3.2 


207538 




















46.6 


± 


6.7 


57.4 ± 5.4 


46.6 


± 


8.3 


28. 0± 3.2 


224151 


23.4 


± 


2.7 








28.9 


± 


2.0 








28.4 ± 1.5 


23.2 


± 


1.5 


5.8 ± 0.7 


303308 


22.6 


± 


2.4 


25.3 ± 2.5 














24.6 ± 1.5 


24.2 


± 


2.8 


5.3 ± 1.1 


308813 


21.3 


± 


2.5 


26.8 


± 3.9 


38.4 


± 


2.4 


35.7 


± 


3.6 


31.6 ± 2.6 


24.4 


± 


2.4 


4.0 ± 1.5 
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Table 5. Absorption Lines and Oscillator Strengths 



HD Lines C I*lines 

Wavelength (A) log(fA) a (cm) Wavelength (A) log(fA) b (cm) 



959.7968 


-6.8520 


945.338 


-5.842 


975.5524 


-6.7190 


1157.4056 


-7.759 


1011.4439 


-6.5769 


1157.7697 


-6.993 


1021.4436 


-6.5867 


1158.5443 


-8.107 


1054.2800 


-6.7632 


1158.6744 


-7.770 


1066.2636 


-6.9114 


1158.7321 


-7.702 


1105.8335 


-8.0849 


1276.7498 


-7.405 






1279.0562 


-7.676 






1279.8907 


-6.737 






1280.4043 


-7.249 






1280.5975 


-7.045 






1656.2672 


-6.004 






1657.3792 


-6.228 






1657.9068 


-6.103 



a HD log(fA) values from Abgrall & Roueff (2006) 
b C I* log(fA) values from Morton (2003) 



Table 6. 


1105.73 A C I* f- value 


Star 


Wll05. 73 


log N(C I*) 


f-value 




(mA) 


(cm -2 ) 


(io- 3 ) 


HD 12323 


16 ± 3 


14.52 


6.011;! 


o Per 


17 ± 4 


14.88 




HD 207538 


32 ± 5 


14.93 


9-o±i:| 


HD 210839 


29 ± 10 


14.85 




Weighted Avg 






6.2t?J 
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Table 7. HD Column Density 



Star (HD) 


Single Component COG 


Multiple Component COG 




log N(HD) 6-value 


log N(HD) 



24534 


15.88 


± 


0.40 a 


1.8 


± 


0.2 


15.71 


± 0.07 


27778 


15.89 


± 


1.10 


2.2 


± 


0.6 


15.84 


± 0.12 


53367 


15.74 


± 


0.34 


4.0 


± 


0.7 


16.23 


± 0.17 


73882 


15.81 


± 


1.18 


2.9 


± 


1.1 


16.03 


± 0.34 


110432 


15.43 


± 


0.12 a 


2.3 


± 


0.2 


15.30 


± 0.08 


147888 


15.39 


± 


1.47 


1.0 


± 


0.1 


15.30 


± 0.15 


149404 


15.74 


± 


0.12 a 


1.1 


± 


0.3 


15.50 


± 0.18 


185418 


15.51 


± 


0.32 a 


3.7 


± 


0.6 


15.87 


± 0.32 


192639 


15.57 


± 


0.22 a b 


4.0 


± 


0.4 


15.10 


± 0.14 


199579 


15.08 


± 


0.09 a 


6.5 


± 


0.7 


15.09 


± 0.02 


206267 


15.72 


± 


0.14 a 


3.4 


± 


0.2 


15.54 


± 0.10 


207198 


16.12 


± 


0.47 a 


3.4 


± 


0.1 


15.84 


± 0.08 


207538 


15.83 


± 


0.25 a 


3.6 


± 


0.5 


15.82 


± 0.09 



12323 


15. 


.02 


± 


0. 


18 


5. 


.4 ± 0. 


.8 


15558 


15. 


.64 


± 


0. 


32 


6. 


.4 ± 1. 


1 


45314 


15. 


.30 


± 


0. 


37 a 


4. 


.9 ± 0. 


7 


74920 


15. 


.33 


± 


0. 


04 a 


3 


.2 ± 0. 


1 


90087 


14. 


.30 


± 


0. 


07 


5. 


.7 ± 2. 


6 


91651 


13. 


.43 


± 


0. 


.45 




>0.6 




91824 


14. 


.17 


± 


0. 


07 a 




>4.8 




93204 


14. 


.22 


± 


0. 


10 a 




>4.2 




93205 


14. 


.28 


± 


0. 


07 




>5.7 




93206 


14. 


.19 


± 


0. 


10 


4. 


.1 ± 2. 


.7 


93222 


14. 


.19 


± 


0. 


12 




>3.2 




94493 


14. 


.24 


± 


0. 


32 




>1.3 




101131 


14. 


.85 


± 


0. 


19 a 


2 


.4 ± 0. 


3 


101190 


15. 


.10 


± 


0. 


12 


2 


.2 ± 0. 


1 


101413 


15. 


.01 


± 


0. 


31 


2. 


.0 ± 0. 


3 


101436 


15. 


.36 


± 


0. 


17 


1 


.4 ± 0. 


5 


104705 


11 


.26 


± 


0. 


53 


1. 


.5 ± 1. 


2 


116852 


14. 


.15 


± 


0. 


13 


2 


.0 ± 1. 


5 


148422 


15. 


.07 


± 


0. 


28 


2 


.0 ± 0. 


3 


152233 


11 


.95 


± 


0. 


19 


2 


.9 ± 0. 


3 


152248 


15. 


.15 


± 


0. 


18 a 


2 


.6 ± 0. 


3 


152723 


15. 


.01 


± 


0. 


13 a 


4. 


.5 ± 0. 


5 


161807 


11 


.12 


± 


0. 


07 a 




>4.5 




177989 


11 


.92 


± 


0. 


22 a 


2 


.6 ± 0. 


.4 


201345 


14. 


.00 


± 


0. 


04 




>5.8 




224151 


11 


.91 


± 


0. 


09 a 


2 


.4 ± 0. 


2 


303308 


11 


.90 


± 


0. 


11 


2 


.2 ± 0. 


3 


308813 


14. 


,18 


± 


0. 


14 


5. 


.2 ± 2. 


1 



a Assumed a good fit then scaled errors to a reduced chi square of one. 

b The single COG is a better solution and agrees with the Profile fitting 
column densities within errors while the multiple component COG is not 
within errors and misfit the data (see Figure [3). 
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Table 8. Column Density Comparison 



Star 


Profile Fitting 


Curve of Growth 




log N(HD) 


b- value 


log N(HD) 


b- value 




(cm- 2 ) 




(cm- 2 ) 




90087 


14.36 ± 0.05 


3.59 


14.30 ± 0.07 


5.7 


101436 


15.12 ± 0.08 


1.33 


15.36 ± 0.17 


1.4 


177989 


14.92 ± 0.07 


2.59 


14.92 ± 0.22 


2.6 


110432 


15.32 ± 0.04 


2.32 


15.43 ± 0.12 


2.3 


192639 


15.52 ± 0.07 


3.70 


15.57 ± 0.22 


4.0 



Table 9. Correlation Coefficients 



Quantities 




P 


# 


Slope 


log N(HD) vs. log N(H 2 ) 


0.94 


0.89 


0.96 


55 


1.25 ± 0.03 


log N(HD) vs. log N(CO) 


0.93 


0.77 


0.98 


13 


0.49 ± 0.02 


log N(HD) vs. logN(H Tot ) 


0.67 


0.48 


0.80 


48 


1.80 ± 0.05 


log N(HD) vs. E(B-V) 


0.78 


0.64 


0.87 


51 


3.94 ± 0.09 


log N(H 2 ) vs. E(B-V) 


0.75 


0.59 


0.85 


51 


2.34 ± 0.04 


log N(HD) vs. Ay 


0.56 


0.27 


0.75 


35 




log N(H 2 ) VS. Ay 


0.53 


0.23 


0.73 


35 




log N(Fe II/H Tot ) vs. log N(HD) 


-0.78 


-0.92 


-0.44 


15 


-0.30 ± 0.07 


log N(Fc II/H Tot ) vs. log N(HD/H Tot ) 


-0.72 


-0.90 


-0.33 


15 


-0.33 ± 0.09 


log N(Fe II/H Tot ) vs. log N(HD/2H 2 ) 


-0.80 


-0.93 


-0.48 


15 


-0.85 ± 0.17 


Correlations not shown 


log N(HD) vs. log N(H I) 


0.34 


0.07 


0.57 


48 




log N(H 2 ) vs. log N(H I) 


0.33 


0.05 


0.56 


48 




log N(HD) vs. Rv 


-0.25 b 




-0.54 


35 




log N(H 2 ) vs. Rv 


-0.26 b 




-0.54 


35 




f H 2 vs. log N(HD) 


0.74 


0.58 


0.85 


48 




Shi vs. log N(H 2 ) 


0.81 


0.68 


0.89 


48 





a r is the Pearson correlation coefficient, p is the 95% confidence interval (uncertainties 
in the correlation coefficients), # is the number of data points used in the correlation, 
and Slope is the slope of the weighted least-squares line fit. 

b Correlations not 95% significant 
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Fig. I. — Carbon I* Curve of growth. 
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Fig. 2. — Single Component (left) and Multiple Component (right) Curves of Growth for 
HD 24534 through HD 149404. 
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Fig. 3. — Single Component (left) and Multiple Component (right) Curves of Growth for 
HD 185418 through HD 207538. 
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Fig. 4. — Column Density solutions from a multiple component curve of growth compared 
to a single component curve of growth. 



-26- 



-4.0 
-4.5 
-5.0 

-5.5 

-6.0 

-4.0 

-4.5 

-5.0 

-5.5 
-6.0 

-4.0 
-4.5 
-5.0 

-5.5 

-6.0 

-4.0 

-4.5 

-5.0 

-5.5 
-6.C 

-4.0 

-4.5 

-5.0 

-5.5 
-6.0 

-4.0 
-4.5 
-5.0 

-5.5 

-6.0 

-4.0 

-4.5 

-5.0 

-5.5 
-6.0 



HD 12323 


HD 1 5558 


HD 45314 


log N(HD) = 15.02 
/ b = 5.4 km s" 


/ log jV(HD) = '5.64 
/ b = 6.4 km s _1 


/ log iV(HD) = 15.30 
" / b = 4.9 km s"' 


u n 7 a Ci o n 

HU /4yzU , 




HU a 1 DO ^ 


/ log TV(HD) = 15.33 
" / b = 3.2 km s"' 


/ log yV(HD) = 14.30 
/ b = 5.7 km s~' 


oJf log iV(HD) = 13.43 
* b a 0.6 km s"' 


HD 9 1824 

/ log iV(HD) = 14.17 
~ / b > 4.8 km s~' 


^ HD 93204 

/ log yV(HD) = '4.22 
/ b s 4.2 km s"' 


^ HD 93205 

/ log iV(HD) = 14.28 
V b s 5.7 km s _1 


^ HD 93206 


_ HD 93222 


_ HD 94493 


/ log iV(HD) = 14.19 

" / b = 4.1 km s" 1 


/ log yV(HD) = 14.19 
/ b a 3.2 km s"' 


log JV(HD) = 14.24 
/ b s 1 .3 km s~' 


HD 101 131 


HD 101190 


HD 101413 


s/ log iV(HD) = 14.85 
b = 2.4 km s"' 


log yV(HD) = 15.10 
/ b = 2.2 km s~' 


>R log yV(HD) = 15.01 
/ b = 2.0 km s"' 


HD 101436 


HD 104705 


_ HD 1 16852 


bg iV(HD) = 15.36 
/ b = 1.4 km s"' 


log iV(HD) = '4.26 
/ b = 1 .5 km s _1 


/' log iV(HD) = 14.15 
"/ b = 2.0 km s"' 


^ HD 148422 


^ HD 152233 


HD 152248 


X log /V(HD) = 15.07 
/ b = 2.0 km s" 1 


>K log yV(HD) = '4.95 
/ b = 2.9 km s _1 


log iV(HD) = 15.15 
/ b = 2.6 km s"' 



9 10 



7 8 9 10 11 6 

NfA) (NfA in cm" 1 



10 11 



Fig. 5. — Single Component Curves of Growth with column densities and 6-values for HD 
12323 through HD 152248. 
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Fig. 6. — Single Component Curves of Growth with column densities and b- values for HD 
152723 through HD 308813. 




Fig. 7. — Profile fit of HD 110432's seven absorption lines, simultaneously solving for column 
density (log N(HD) = 15.32) and 6-value (2.32). The velocity is set to the rest wavelength 
due to varying velocity solutions from different detector segments. 
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Fig. 8. — Correlation plots of log N(HD) versus log N(H 2 ). The line is a weighted least- 
squares fit. Square plot symbols represent Copernicus data analyzed by Lacour et al. 2005, 
diamonds represent this study, asterisk and upper limit represent Copernicus data reported 
by Spitzer et al. 1974. 
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Fig. 9.— Correlation of log N(CO) with log N(HD). CO data from Burgh et al. 2007 
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Fig. 10. — Correlation plots of log N(HD) versus log N(total Hydrogen). Square plot symbols 
represent Copernicus data analyzed by Lacour et al. 2005, diamonds represent this study. 
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Fig. 11. — log N(HD) and log N(H 2 ) versus reddening, E(B-V). Symbols same as Figure fTUl 
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Fig. 12. — log N(HD) and log N(H 2 ), versus extinction (Ay). Symbols same as Figure [TU1 
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Fig. 13. — Correlations with iron depletions, a density indicator. Iron data in Table [T] 
originally reported by Jensen Sz Snow 2007. 



